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Abstract. We have verified the absolute timing capabilities of the high-energy instruments aboard INTEGRAL, 
i.e. the imager IBIS, the spectrometer SPI and the X-ray monitor JEM-X. Calibration observations of the Crab, 
contemporaneous with the Rossi X-ray Timing Explorer (RXTE), have been used to measure the absolute phase 
of the main pulse of the Crab pulse profile using the same Jodrell Bank radio ephemeris. The three INTEGRAL 
instruments and RXTE give within the statistical and systematic uncertainties consistent results: The X-ray main 
pulse is leading the radio pulse by 280 ± 40/is. Also the shapes of the X-ray pulse profiles as measured by the 
different instruments are fully consistent with each other. In addition, we present the first measurement of the 
absolute phase of the main pulse at 7-ray energies above 30 MeV using data from the EGRET instrument aboard 
the Compton Gamma- Ray Observatory: The 7-ray main pulse is leading the radio one by 241 ± 29/is, consistent 
with the value for the X-ray main pulse. Comparing absolute arrival times at multiple frequencies gives important 
constraints to models explaining the production of non-thermal emission in magnetospheres of rotation powered 
neutron stars. 
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1. Introduction 

The high-energy instruments aboard INTEGRAL 
(Winkler et al. l2003[) are designed to allow timing studies 
of e.g. radio pulsars with periods down to the very short 
periods of millisecond pulsars. This requires very precise 
and stable timing information, particularly for studies at 
the high-energy end of the INTEGRAL window where 
the counting rates are low and long integration times 
required. An important aspect in timing studies is the 
measurement of absolute arrival times. Particularly in 
multiwavelengths studies of pulsars, the comparison of 
absolute arrival times is an important diagnostic in the 
detailed modelling of the structures of the production 
sites in magnetospheres of rotation-powered neutron 
stars. 



We used INTEGRAL calibration observations of the 
Crab pulsar (PSR B0531+21) and contemporaneous ob- 
servations with the Proportional Counter Array (PCA, 
Jahoda 1994) on the Rossi X-ray Timing Explorer 
(RXTE) to verify the correct absolute timing calibration 
for the imager IBIS (Ubertini et al. I2003[) . the spectrom- 
eter SPI (Vedrenne et al. 2005)) and the X-ray monitor 
JEM-X (Lund et al. [2003 ) . This was done by measuring 
the absolute phase of the main pulse of the Crab pulse pro- 
file with the INTEGRAL instruments and the PCA. All 
data have been analyzed with the same Jodrell Bank radio 
ephemeris, reducing systematic effects due to e.g. disper- 
sion measure variations causing systematic variations in 
the measured arrival times. 



Send offprint requests to: L. Kuiper, 
e-mail: L.M.Kuiper@sron.nl 

* Based on observations with INTEGRAL, an ESA project 
with instruments and science data centre funded by ESA 
member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain), Czech Republic and 
Poland, and with the participation of Russia and the USA. 



Finally, the arrival time of the main pulse of the Crab 
pulsar has been studied from radio, IR, optical up to X- 
ray wavelengths (e.g. Ray et al. I2003|) . We present here a 
first estimate of the 7-ray (30 MeV - 30 GeV) radio delay 
using data from EGRET (Thompson et al. I1993f) aboard 
the Compton Gamma- Ray Observatory (CGRO). 
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2. INTEGRAL timing resolution 

In order to stay within the limited telemetry budget of 
INTEGRAL, the event timing information measured by 
the instruments has been degraded on board. The accu- 
racy of the On Board Time (OBT) for IBIS/ISG RI (op- 
timized for energies ~ 15 - 200 keV, Lebrun et al. I2003[) 
is thus 61/xs, for IBIS/PICsIT (optimized for energies 175 
keV - 10 MeV, Di Cocco et al. |2~UU3|) is 64/zs (for IBIS 
the least significant 8 bits are cleared), for SPI 102/j,s (the 
least significant 11 bits are cleared) and for JEM-X 122/is 
(the readings of the Digital Front End Electronics clock). 
There are additional smaller uncertainties affecting the 
event timing accuracy, e.g. the OBT might have an uncer- 
tainty up to about 16/is, uncertainties in the orbit predic- 
tions can cause similar uncertainties and for pulsar stud- 
ies the conversion to the arrival time at the Solar System 
Barycentre (SSB) will again add a similar error. These un- 
certainties affect the timing of single events and cause a 
broadening of derived timing signatures. We estimate that 
the resulting timing resolutions for pulsar timing studies 
using the arrival times at the SSB are 90/is for IBIS, 130/is 
for SPI and 150/xs for JEM-X (3tr accuracy). As a first 
check, we used the Crab calibrations to compare the pulse 
profiles measured by the INTEGRAL instruments with 
the RXTE profile. A more accurate calibration should be 
performed in the future by studying the profiles of mil- 
lisecond pulsars. 

3. Observations and data 

INTEGRAL observed the Crab from 7 February 2003 
(revolution 39) up to and including 27 February 2003 (rev- 
olution 45). The spacecraft and instruments were in dif- 
ferent modes during this calibration period. For our in- 
vestigations we need photon-by-photon modes in which 
each event is time tagged. The IBIS/ISGRI detector, SPI 
and JEM-X were all in photon-by-photon mode during 
revolution 42, roughly in the middle of the Crab obser- 
vations. We used for these three instruments the data 
from revolution 42, collected at the REDU groundsta- 
tion, which provided ample statistics. The IBIS/PICsIT 
detector, unfortunately, cannot routinely be configured in 
photon-by-photon mode due to PICsIT's high count rate 
and the tight telemetry budget of INTEGRAL. In revolu- 
tions 39 and 40, PICsIT was part of the time commanded 
in photon-by-photon mode allowing a test of the timing 
capabilities. We used these data collected at REDU. 

The RXTE PCA (Jahoda ITU9ljl routinely monitors 
the Crab pulsar, and had one observation on 13 February 
2003, roughly in the middle of our Crab calibrations, and 
one on 27 February 2003, at the end of the INTEGRAL 
observations. We retrieved these data from the High 
Energy Astrophysics Science Archive Research Center 
(HEASARC) at NASA/Goddard Space Flight Center. 

EGRET, sensitive to high-energy 7-rays from 30 MeV 
to 30 GeV (Thompson et al. H533)l . had the Crab 20 times 
in its field-of-view with high signal-to- noise over a ~ 7.5 



Table 1. Radio ephemcris for PSR B0531+21 (monthly 
monitoring at Jodrell Bank) valid during the INTEGRAL 
Crab calibration. 



Parameter Value 

Val. range (MJD) 52 671-52 699 

to (MJD; TDB) 52 685.0 

u (Hz) 29.8092705122612 

i> (10~ 10 Hz s" 1 ) -3.73661 

i>o (10~ 20 Hz s~ 2 ) 2.01 

RMS of ToA fit 23.5/ns 

Phase of radio MP at zero epoch (to) 0.2838 

a 2 ooo 05 h 34 m 31 s .972 

5 20 oo 22°00'52".07 



year baseline starting in April 1991 (details of the EGRET 
observations of the Crab are given in Kuiper et al. (.2001) . 
We used these data to study the 7-ray radio delay of the 
main pulse. 

4. Timing analysis 

The Jodrell Bank radio telescope (UK) performs monthly 
monitoring of PSR B0531+21, providing us with a Crab 
radio ephemcris valid for INTEGRAL revolutions 39- 45, 
including the dates of the RXTE observations (see Table 
I)- 

4.1. RXTE analysis and results 

In the timing analysis of the RXTE PCA data of 13 and 
27 February 2003 the barycentered event arrival times 
were folded with the ephemeris in Table 1, providing us 
with accurate X-ray pulse profiles in absolute phase (see 
Fig. ^). The PCA events are time tagged with a 1/xs ac- 
curacy with respect to the spacecraft clock and with an 
absolute time accuracy of 5 - 8/is with respect to UTC. 
The maximum of the main radio pulse is at phase 0.0, 
the fiducial point in the radio template. For both data 
sets (13 and 27 February), we determined the absolute ar- 
rival time/phase of the main X-ray pulse by fitting the 
profile with an asymmetric Lorentzian plus a constant 
background. The fits were perfect with the X-ray maxima 
leading the radio pulse by 248 ± 7 [is and 289 ± 7^s, respec- 
tively. The errors are ler, statistical only. The difference 
between the two estimates is a measure of the systematic 
uncertainty in the radio parameters, e.g. due to dispersion 
measure variations, and is consistent with earlier RXTE 
estimates of 302 ± 67/is (Rots et al. 120001 and Tennant et 
al. l200U and also with a recent estimate of ~ 300/xs using 
XMM-Newton data (Kirsch et al. EHOHl- 

4.2. INTEGRAL analysis and results 

For the analysis of the INTEGRAL data, we selected all 
time-tagged events received via the REDU groundstation 



L. Kuiper et al.: Absolute timing with IBIS, SPI and JEM-X aboard INTEGRAL 



1.22xl0- a 



7.76X10- 3 



3.38X10- 3 



-l.OlxlO- 3 I 

0.7 



1.22x10 



1.1 1.3 

Phase 




7.76 



3.36 



-1.01x10 




Phase 



1.22xl0- a 



7.76xl0- 3 



3.38x10-3 



1.22xl0- a 



7.06x10-3 - 



1.97x10-3 - 



-3.12x10- 





0.7 



0.9 



1.1 



1.3 



1.5 



1.7 



Phase 



Fig. 1. Comparison in absolute phase of Crab pulse profiles measured by INTEGRAL and RXTE: Upper-left JEM-X 
(6 - 23 keV, orange); upper-right IBIS/ISGRI (20 - 50 keV, red); lower-left SPI (20 - 50 keV, green); lower-right 
IBIS/PICsIT (180 - 580 keV, blue). The RXTE profile (2 - 30 keV, black) is produced with the same ephemeris and 
shown in each figure. 



during revolution 42. For both JEM-X monitors Pulse 
Height Analyzer channels 128 - 192 were selected, corre- 
sponding at the time of the Crab calibration for JEM-X1 
to energies 9-23 keV, and JEM-X2 6-15 keV (a com- 
plementary timing study of the Crab pulsar with JEM-X 
is reported by Brandt et al. IM^ . For IBIS/ISGRI and 
SPI events were used with energies 20 - 50 keV. Finally, 
for IBIS/PICsIT we had ~ 214 ks exposure in photon-by- 
photon mode in revolutions 39 and 40, and selected events 
with energies 180 - 580 keV (PHA 25 - 75). Unfortunately, 
the latter data set has low statistics. 

The arrival times of the events at the spacecraft are 
registered in On Board Time (OBT), which is converted 



on ground to Terrestrial Time, expressed in INTEGRAL 
Julian Date (IJD) (see Walter et al. ED® - Th e correla- 
tion between the OBT and IJD needs to take into account 
delays on board for each of the instruments, delays within 
the ground stations and due to light travel time. A full 
description with correction values is given in Walter et 
al. (2003). In this work we applied all these corrections 
valid during the Crab calibrations to obtain absolute time: 
Offset added to UTC 964^s; delays ESA REDU groun- 
station — 103/xs and NASA DSS-16 groundstation 42/is; 
delays per instrument derived from ground calibrations, 
IBIS 111/zs, SPI 134/iis and JEM-X 185/zs. 
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The arrival times of the INTEGRAL time-tagged 
events were also converted to the arrival time at the Solar 
System Barycentre and folded with the Crab ephemeris in 
Table 1. The resulting Crab profiles in absolute phase (i.e. 
the main radio pulse is at phase 0.0) are shown in Fig. ^ 
with the profiles of JEM-X1 and JEMX-2 summed in one 
figure. The excellent agreement with the RXTE PCA pro- 
files is evident for the JEM-X data. At higher energies the 
INTEGRAL data display the known evolution of the Crab 
pulse profile with increasing energy: 1) a higher level of 
bridge emission between the peaks for IBIS and SPI com- 
pared to RXTE and JEM-X, 2) the relative intensity of 
the first pulse is higher for IBIS and SPI than for RXTE 
and JEM-X (see e.g. Kuiper et al. I2001[) . On this scale, 
we do not see any broadening of the INTEGRAL pro- 
files relative to the RXTE profile. For comparison, the 
half-width-half-maximum of the first Crab pulse is about 
600^s. 

In order the compare the absolute times more accu- 
rately, we determined for JEM-X1, JEM-X2, IBIS/ISGRI 
and SPI the arrival phases of the main (first) pulses, simi- 
lar to what was done for RXTE. In each case the profile of 
the main pulse was fitted with an asymmetric Lorentzian 
plus constant background and the phase of the maximum 
was determined. The uncertainties in the latter have a 
small statistical error, but there is also a systematic un- 
certainty which depends on the choice of the fit func- 
tion. In order to get an estimate of the magnitude of 
the systematic uncertainties, we also made Fourier fits 
with 40 harmonics. Fig. |3 shows as an example the fits 
for IBIS/ISGRI. The phase difference between the two fit 
maxima is 15/iS, representative for the systematic uncer- 
tainties for all instruments. 

For the Fourier fits, we find the following X-ray ra- 
dio delays: JEM-X1 284 ± 17/xs, JEM-X2 300 ± 14/zs, 
IBIS/ISGRI 285±12^s, SPI 265±23/is. The errors are la, 
statistical only. Taking the systematic uncertainties into 
account, the average value for the measured X-ray radio 
delay is ~ 280 ± 40/iS. It can be noted that, firstly, the 
INTEGRAL instruments give fully consistent results, and, 
secondly, the average INTEGRAL value is fully consistent 
with the RXTE PCA values derived using the same radio 
ephemeris, even though the two X-ray measurements em- 
ploy very different means of measuring time and satellite 
position and thus have different systematic error budgets. 

Due to the tight telemetry budget, it is foreseen 
that IBIS/PICsIT will hardly be configured in photon- 
by-photon mode. For timing studies observers can select 
the spectral-timing mode with highest time resolution of 
0.97656 ms. To verify the correct performance, PICsIT 
was in revolutions 40 and 41 also in spectral-timing mode 
for 345 ks. We produced the Crab pulse profile shown in 
Fig. |21 which has a significance of 16o\ The profile is con- 
sistent with those shown in Fig. ^ The intense bridge 
emission between the peaks is expected for the selected 
high energy range (band 3, 260 - 364 keV). 
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Fig. 2. Main pulse of the Crab profile measured by 
IBIS/ISGRI. Asymmetric-Lorentzian fit (green) and 
Fourier fit with 40 harmonics (red). Phase shift between 
the maxima of the two fits is 15/^s. The x-axis is not in 
absolute phase, the known corrections factors (see text) 
are not yet applied. 
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Fig. 3. Crab pulse profile in absolute phase measured by 
IBIS/PICsIT in spectral-timing mode for band-3 (260 - 
364 keV) using 345.2 ks exposure (blue histogramme in 
30 bins). The RXTE profile (black, 2-30 keV) is shown 
for comparison. 

4.3. Arrival time Crab main pulse at 7-ray energies 
above 30 MeV 

For the 20 EGRET observations with high-statistics Crab 
data, we constructed for all observations the high-energy 
(30 MeV - 30 GeV) pulse profiles, following Kuiper et al. 
2001 It should be noted that different ephemerides had 
to be used over the CGRO mission lifetime, such that the 
different profiles can not be directly summed in absolute 
phase. Due to e.g. variations in the dispersion measure, de- 
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Fig. 4. Arrival phase of the main pulse vs time for photon 
energies 30 MeV - 30 GeV derived from CGRO/EGRET 
data. The Ha uncertainty interval assuming a constant 
arrival phase is indicated by the hatched region. The high- 
energy 7-rays arrive 241 ± 29/Lts earlier than the fiducial 
point in the radio template. 

rived arrival phases will scatter around the genuine value. 
The main peak of the high-energy profile is very symmet- 
ric, because the bridge emission is very weak at these high 
energies. Therefore, we determined for each observation 
the arrival phase of the main pulse by fitting a symmetric 
Lorentzian profile plus a constant background to the mea- 
sured event distribution in a relatively broad (0.3) phase 
window centered around phase 0. In Fig. 0] these arrival 
phases are shown as a function of time over a ~ 7.5 year 
baseline. The error bars on the individual datapoints are 
composed of a statistical contribution from the profile fit 
and of a systematic component of about 100/zs (lex) re- 
flecting the uncertainties in the dispersion measure from 
the radio pulsar timing observations (see explanatory sup- 
plement README.txt in the Princeton CGRO database 
directory). The arrival phases have been fitted by a con- 
stant value yielding as optimum value — 0.0072±0.0009 for 
an acceptable xt value of 1.59 for 19 degrees of freedom. 
This offset in the phase domain translates to an offset 
of 241 ± 29/is in the time domain and indicates that the 
high energy 7-ray photons from the main pulse also ar- 
rive earlier than the maximum of the main radio pulse. 
This value of 241/zs is significantly larger than the quoted 
CGRO absolute timing uncertainty of better than 100/is 
and is fully consistent with the results obtained above from 
the INTEGRAL and RXTE X-ray data. 

5. Discussion 

We have determined with all high-energy instruments 
aboard INTEGRAL, i.e. the imager IBIS, the spectrom- 
eter SPI and the X-ray monitor JEM-X, the X-ray radio 
delay between the absolute arrival times of the X-ray and 
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Fig. 5. Main pulse of the Crab profile in absolute phase as 
measured by IBIS/ISGRI (20 - 50 keV; full line) compared 
to the Greenbank 812.5 MHz radio profile (broken line; 
Lundgren 1994). 

radio main pulses in the profile of the Crab pulsar, PSR 
B0531+21. Within the statistical accuracies these instru- 
ments measured the same delay, on average ~ 280 ± 40/is 
(the error includes estimates of statistical and systematic 
uncertainties), the X-ray pulse leading the radio one. This 
value is consistent with the delays we measured twice with 
the RXTE/PCA during the period of the INTEGRAL 
Crab calibrations, namely 248 ± 7/is and 289 ± 7/xs, sta- 
tistical errors only. These results show that INTEGRAL 
meets its specifications and can measure absolute arrival 
times with accuracies down to about 40/j.s. Fig. [5] shows 
the X-ray main pulse measured by ISGRI for energies 20 - 
50 keV together with its radio counterpart at 812.5 MHz 
measured with the Greenbank 140 ft telescope (Lundgren 
1994). In addition to the phase shift, it is evident that the 
profile shapes are very different. The narrow radio profile 
is shifted but contained in the phase range covered by the 
X-ray profile. 

We also measured the delay of the arrival time of the 
main pulse in the high-energy (30 MeV - 30 GeV) 7-ray 
profile relative to the main radio pulse, using 20 Crab ob- 
servations by CGRO/EGRET over a ~ 7.5 year baseline: 
241 ± 29/is. This value appears to be consistent with the 
delay measured for the main X-ray pulse and is an ad- 
ditional constraint in the three-dimensional modelling of 
the production processes and sites in the magnetospheres 
of rotation powered pulsars. 

In Kuiper et al. 120011 presenting for the Crab a co- 
herent high-energy picture from soft X-rays up to high- 
energy 7-rays, it is argued that models with an outer gap 
are more successful in reproducing the overall character- 
istics of PSR B0531+21 as measured from the radio up to 
the 7-ray domain, than polar cap models. These are the 
two competing model scenario's attempting to explain the 
high-energy radiation from magnetized rotating neutron 



6 



L. Kuiper et al.: Absolute timing with IBIS, SPI and JEM-X aboard INTEGRAL 



stars. The known large magnetic inclination angle a, i.e. 
the angle between the pulsar spin axis and magnetic mo- 
ment axis, and aspect angle £, i.e. the angle between the 
spin axis and the observer's line of sight, which are for 
the Crab both estimated from radio and optical observa- 
tions to be ~ 60°, cause for polar cap scenarios problems 
to reproduce the structure of the Crab light curve (see 
e.g. discussion in Kuiper et al. I2001|) . Fig. shows the 
phase aligned total Crab profiles for the energy windows 
we are addressing in this work (for a detailed presenta- 
tion of the frequency-dependent behaviour of the average 
radio pulse profile, see Moffett & Hankins 1996). Outer 
gap scenarios can for these large geometrical angles a and 
£ account for the measured double peaked profile with 
bridge emission in the X-ray and soft 7-ray band (e.g. 
Chiang & Romani 1994), and are also succesful in repro- 
ducing characteristics as the observed polarization posi- 
tion angle swings of the Crab pulsar at optical wavelengths 
(Romani & Yadigaroglu 1995). The radio precursor, vis- 
ible in the 430 MHz profile in Fig. H3 a, originates in the 
outer gap models from near the polar cap. The phase shift 
between the high-energy main pulse (X-rays and 7-rays) 
and the main radio pulse is an additional detail to be ac- 
counted for. More recently, Cheng, Ruderman & Zhang 
2000. extended the outer magnetospheric gap model to a 
detailed three dimensional model with photon emission 
moving outward and inward from regions in outer gaps 
above both poles, the gaps being limited along the az- 
imuthal direction by e^- pair production of inward-flowing 
photons from the outer gaps. In their model the main X- 
ray/7-ray pulse is produced high in the magnetosphere at 
0.8 < r/Ri c < 1.0 with Ri c the light cylinder radius. In 
order to account for the delay of the radio main pulse, the 
production site of radio emission in the narrow outer gap 
should be slightly shifted inwards to the neutron star with 
respect to the production site for the high-energy emiss- 
sion, but still contained in the latter. It is expected that 
the intensity weighting along the gap surface differs for dif- 
ferent energy bands, but the non-thermal radio emission is 
likely produced in a narrower region of the thin outer gap 
along the last closed field line. More detailed modelling 
should show whether this can explain the profiles shown 
in Fig. \5\ 

Acknowledgements. We wish to thank the High Energy As- 
trophysics Science Archive Research Center (HEASARC) at 
NASA/Goddard Space Flight Center for maintaining its online 
archive service which provided the RXTE data used in this re- 
search. LF acknowledges financial support by the Italian Space 
Agency (ASI) and the hospitality of the ISDC during part of 
this work. 

References 

Brandt, S., Budtz-j0rgensen, C, Lund, N., et al. 2003, A&A, 
this volume 

Cheng, K.S., Ruderman, M., Zhang, L., 2000, ApJ 537, 964 
Chiang, J., Romani, R.W., 1994, ApJ 436, 754 
Di Cocco, G., Caroli, E., Celesti, E., et al. 2003, A&A, this 
volume 





430 MHz 

LJ 




a 

LJ 




B12.5 MHz 


u 


b 




20-50 keV 


U 




>100 MeV 


a 



0.0 0.5 1.0 1.5 2.0 



Phase 

Fig. 6. Crab profiles in absolute phase: a) radio 430 MHz 
(Arecibo: Moffett & Hankins H^rJjl . b) radio 812.5 MHz 
(Greenbank: Lundgren I1994[> . c) X-rays 20 - 50 keV 
(IBIS/ISGRI), d) 7-rays with energies above 100 MeV 
(EGRET) ). 

Jahoda, K. 1994, AAS Meeting 26, 894 

Kirsch, M.G.F., Becker, W., Benlloch-Garcia, S., et al. 2004, 
to appear in the proceedings of the SPIE 2003 conference 

Kuiper, L., Hermsen, W., Cusumano, G., et al. 2001, A&A, 
378, 918 

Lebrun, F., Leray, J. P., Lavocat, P., et al. 2003, A&A, this 
volume 

Lund, N., Brandt, S., Budtz-J0rgensen, C, et al. 2003, A&A, 
this volume 

Lundgren, S.C. 1994, PhD Thesis, Cornell University 

Moffett, D.A. & Hankins, T.H. 1996, ApJ 468, 779 

Ray, PS., Wood, K.S., Wolff, M.T., et al. 2003, Am. Astrono. 

Soc. HEAD meeting 35, 2213 
Romani, R.W., Yadigaroglu, LA. 1995, ApJ 438, 314 
Rots, A.H., Jahoda, K, Lyne, A.G., 2000, Am. Astron. Soc. 

HEAD meeting, 32, 3308R 
Tennant, A.F., Becker, W., Juda, M., et al. 2001, ApJ 554, 

L173 

Thompson, D.J., Bertsch, D.L., Fichtel, C.E., et al. 1993, ApJS 
86, 629 



L. Kuiper et al.: Absolute timing with IBIS, SPI and JEM-X aboard INTEGRAL 



7 



Ubertini, P., Lebrun, F., Di Cocco, G., et al. 2003, A&A, this 
volume 

Vedrenne, G., Roques, J.-R, Schonfelder, V., et al. 2003, A&A, 
this volume 

Walter, R., Favre, P., Dubath, P., et al. 2003, A&A, this volume 
Winkler, C, Courvoisier, T.J.-L., Di Cocco, G., et al. 2003, 
A&A, this volume 



